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Abstract 

Over the recent past, the galactic environment of the Sun has differed substan- 
tially from today. Sometime within the past ^130,000 years, and possibly as 
recent as '^56,000 years ago, the Sun entered the tenuous tepid partially ion- 
ized interstellar material now flowing past the Sun. Prior to that, the Sun was 
in the low density interior of the Local Bubble. As the Sun entered the local 
ISM flow, we passed briefly through an interface region of some type. The low 
column densities of the cloud now surrounding the solar system indicate that 
heliosphere boundary conditions will vary from opacity considerations alone as 
the Sun moves through the cloud. These variations in the interstellar material 
surrounding the Sun affected the paleoheliosphere. 

1 Introduction 

The boundary conditions of the heliosphere at any given point in its history are 
set by the interstellar cloud that happened to surround the Solar System at that 
time. Our variable galactic environment affects the physical properties of the 
heliosphere, and the fluxes of galactic cosmic rays and interstellar byproducts 
reaching the Earth. The work of Hans Fahr has, from the earliest until now, 
provided a strong foundation for our understanding of the effect of the ISM on 
the heliosphere and on Earth |7I ISI IHI ISIIU^ . 

The properties of both astrospheres and the heliosphere are highly respon- 
sive to the boundary conditions supplied by the interstellar medium and the 
interstellar radiation field [ni[Tniini[7TllSniiniESlEI]- The space motion of 
the Sun, when compared to the column densities of interstellar material (ISM) 
towards stars near the Sun, indicates that sometime in the late Quaternary the 
Sun, which has been moving through the very low density region known as the 
Local Bubble, encountered the cluster of local interstellar clouds (CLIC) flowing 
away from the direction of the Scorpius-Centaurus Association [T^llT^imi^ . 
Mediating the interaction between the very low density Local Bubble and the 
tepid CLIC will be a thin interface of some type ISSIIST] . 

Regardless of the Local Bubble plasma pressure, the large contrast (> 10'^) 
between neutral ISM densities in the Local Bubble versus the CLIC must have 
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generated significant changes in the heliosphere, because pickup ions and anoma- 
lous cosmic rays are processed interstellar neutrals. The well-known anticorrela- 
tion between galactic cosmic ray fluxes at Earth and the solar magnetic activity 
cycle, mediated by the heliosphere, suggests that the transition between the 
Local Bubble cavity and CLIC altered both the heliosphere and the galactic 
cosmic ray flux at Earth, with possible implications for the terrestrial climate. 
Heliosphere models have now been constructed for a range of interstellar bound- 
ary conditions jl^l H^l 1531 1^ . The models show that the ISM surrounding the 
Sun during the geological past (or the "palcoLISM") did affect the heliosphere 
during the geological past (or the "paleoheliosphere" ) , and therefore also the 
galactic cosmic ray flux at Earth. 

In this paper, we postulate a geometric model for the CLIC, supplemented 
by equilibrium models that provide the cloud density. The entry epoch of the 
Sun into the CLIC (also known as the Local Fluff) depends on the relative 
Sun-CLIC velocities and the distribution of the CLIC gas. To estimate this 
transition epoch, the distribution of CLIC gas is derived from a simple model 
of the cloud morphology using data for 11° and D° towards nearby stars f i)2.1|l . 
combined with the CLIC density found from photoionization equilibrium models 
of nearby ISM ( W2.2^ . The entry of the Sun into the CLIC is then found from the 
cloudlet velocities in the downwind direction f ti2.1|) . Prior to encountering the 
CLIC, the Sun was in the very low density gas of the Local Bubble interior for 
several million years (331. The limitations of the simple assumptions underlying 
these estimates are discussed briefly in 21 

2 Sun Passage into Very Local ISM 

The epoch when the Sun first encountered the CLIC gas can be estimated from 
data on ISM column densities and velocities towards nearby stars, combined 
with theoretical models that provide the average cloud density. These photoion- 
ization models provide a second important result about the paleoheliosphere, by 
showing clearly that the boundary conditions of the heliosphere, in particular 
the ionization level of hydrogen and the electron density, vary from radiative 
transfer effects alone as the Sun traverses low opacity ISM f i)2.2l (57i i58) . SF02, 
SF06). 

2.1 Local ISM Data and Distribution 

The Sun has recently entered the CLIC, which itself is inhomogeneous. To 
estimate the entry date, we map the CLIC gas by assuming that the distances to 
the CLIC edges are given by iV(II°)/(7T,(II°)), where A^(II°) are column densites 
towards nearby stars. The average space density (n(II°)) is provided by radiative 
transfer models (' ii2.2|l . 

CLIC dynamics indicate an origin related to a superbubble caused by star 
evolution in the Scorpius-Ophiuchus Association JS]. The bulk flow velocity 
vector of the CLIC gas past the Sun is ~-28±4.6 km s~^, from the direction 
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£,6~12°,12° in heliocentric coordinates This vector is based on absorption 
components towards ~60 stars, obtained at resolutions of 0.3-3.0 km s~^. Most 
of the Hyades stars were excluded from the star sample underlying this vector, 
because ISM with a poorly defined relationship to the CLIC is found inside of 
this cluster The LIC is the best understood member of the CLIC, with a 
precisely known heliocentric velocity based on Ulysses He" data of —26.3 ± 0.4 
km s-i from £,b= 3.3°, 15.9° jBS]. The exact value for the CLIC bulk flow 
depends on the underlying star sample, since the flow gradually decelerates 
towards the downwind direction (FS06). Examples of the deceleration are the 
LIC, found in the downwind direction at +26 km s~^, and the Apex cloud, 
found at -35 km s~^ in the upwind direction (HC). 

The velocity vectors of the CLIC, LIC, and the Sun are listed in Table ^ 
in both HC coordinates and the local standard of rest (LSR) velocity frame, 
for both the Standard and Hipparcos-based LSR frame. The uncertainty in the 
LSR occurs because the age distributions of the star samples underlying these 
two LSR vectors were different (H] 021 ■ 

Data on A^(II°) are drawn from observations of D° and 11° towards nearby 
cool and hot stars |SHlEniE2, white dwarf stars |nilESIEnilllESllS3EiElllIll, 
and nearby stars with observations of interstellar Ca+ f|22ll51l^. and references 
therein). The column densities of all velocity components towards a star are 
summed to obtain the total 7V(H°) towards the CLIC surface. For stars with 
D° data, a conversion factor of N{W)/ N(D°)=1.5xlO-^ is used. The total H° 
column densities to the CLIC surface are in the range A^(H°)= 0.3 x 10^^ to 
~ 10^^ cm-2. 

Optical Ca"*" data are also used, however the ratio iV(Ca+)/iV(H°) is highly 
variable. In cold clouds ~99.7% of the Ca is depleted onto dust grains. Both 
depletion and ionization affect Ca+ in warm clouds such as the CLIC; for ex- 
ample Ca++/Ca+>1 if T > 4,000 K and n(e)<0.13 cm'^ The ratio 
A^(II°)/A^(Ca"'") = 10^^ is used, based on the three absorption components ob- 
served in both Ca+ and H° towards a Aql {d=5 pc, ^,6~48°,-9°, e.g., P^ICT 

HH). 

Looking only at stars within 10 pc, average values of (A^(II°)) ~ 10^^ cm ^ 
and (n(II°)) ^ 0.07 cm^^ are found. Based on temperature and turbulence 
data in Redfield and Linsky (2004), objects within 10 pc show a temperature 
and turbulence range of r=l, 700-12, 600 K and ^ = — 5.5 km s~^, with mean 
values 6, 740 ± 2800 K and ^ = 1.9 ± 1.0 km s^^. The ^ variable represents 
deviations from a Maxwellian velocity distribution for the atoms contributing to 
the absorption line components, and as such is a mock turbulence that includes 
unresolved clouds. 

The range of temperatures inferred for the CLIC, combined with the macro- 
turbulence of ±4.6 km s~^, show that the CLIC is inhomogeneous and that the 
boundary conditions of the heliosphere will vary during the next thousands of 
years as the Sun traverses the CLIC. 

The CLIC is fully contained in the nearest '--^35 pc of space. The distances to 
the edges of the CLIC are shown in Figs, ^andl^] These figures are constructed 
from the data listed above, and assuming (n(II°))=0.17 cm^'^ (discussed in 
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Figure 1: The distances to the edge of the CLIC, projected onto the galactic 
plane, as calculated from iV(H°)/(n(H°)) for stars near the Sun. Distances are 
based on observations of H° (dots), D° (dots), and Ca"*" (crosses) (see text). 
The directions i=0°, 90°, 180°, and 270° are labeled. The arrows directed to 
the left show the CLIC motion through the Local Standard of Rest (LSR) based 
on the standard (solid) and Hipparcos (dashed) solar apex motions, while the 
arrows to the right show the motion of the Sun through the LSR based on the 
Hipparcos and Standard apex values (Table 



tl2.2fl . This ISM distribution is assumed to contain no gaps. 

An extended dilute H II region found towards A Sco [7Dj may contribute to 
the excess cloud length towards the galactic center that is indicated by Ca+ data. 
Highly ionized gas at CLIC velocities is found towards both A Sco (^,6=352°,- 
2°) and HD 149499B (d=37 pc, ^,5=330° ,-7°), where for instance N+/N°=1.9 

m- 

2.2 Boundary Condition Variations from Cloud Opacity 

An important influence on the heliosphere while the Sun traverses the low opac- 
ity CLIC is the change in heliospheric boundary conditions caused by ionization 
variations due to the attenuation of photons with energy E >13.6 eV (the ion- 
ization edge of H°). These variations are shown by our radiative transfer (RT) 
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Figure 2: Same as Fig. ^ except that the projection is on a meridian plane 
perpendicular to the galactic plane. 

models [SZ]'^- The CLIC is partially opaque to H-ionizing photons (A < 912 
A) and nearly transparent to He- ionizing photons (A < 504 A). A cloud op- 
tical depth of T ~1 is achieved for logA^(H°)~ 17.2 cm~^ and ^ 17.7 cm^^, 
respectively, at the H° and He° ionization edges. 

A series of models have been constructed to study these opacity effects (SF02, 
SF06). These models employ the CLOUDY radiative transfer code, which incor- 
porates a wide range of physical processes to model a cloud under the conditions 
of photoionization equilibrium |1(J| . The models are constrained by observations 
of the CLIC towards e CMa |57], observations of pickup ions, anomalous cosmic 
rays, and n(He°) inside of the solar system, and interstellar radiation field data 
and models. The radiation field is based on radiation sources affecting the solar 
environment (see SF02), and extends out to soft X-ray energies f ii3.2l Fig. 01 
summarizes the variations in neutral densities that are obtained for equilibrium 
calculations of low density ISM similar to the LIC (see SF02 and SF06). The 
best of these models (models 2 and 8 in SF02) give n(H°)=0. 18-0.21 cm'^ and 
n(H+)=0.1 cm~^ at the solar location. These model results that best predict 
the observed gas densities, for a study still in progress, give a mean value of 

^By "radiative transfer effects", we mean that the spectrum of the radiation field, including 
both point source and diffuse contributions, is substantially modified for energies E > 13.6 
eV as the radiation penetrates more deeply into the cloud. This is the case for the CLIC gas, 
as shown by the relative opacities of the cloud to H- versus He- ionizing photons. 
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Figure 3: The densities of H° and He° , as predicted by radiative transfer models 
of the CLIC, are shown for a range of equilibrium conditions based on models in 
SF02 and SF06. The symbols give the model parameters: symbol colors indicate 
the magnetic field strength in the cloud interface; symbol fill gives gives the 
assumed HI column density; symbol shape gives the assumed temperature of 
the Local Bubble plasma. The stars are special parameter sets which do not fall 
on the grid of model parameters, but rather are chosen to better match the data. 
The range of current values for n(He°), based on Ulysses and pickup ion data, 
are shown as vertical lines [35]. Although a range of n{ll°) values are consistent 
with 7i(He°), the radiative transfer models providing the best agreement with 
ISM observations both inside of the heliosphere and towards nearby stars give 
n(H°)= 0.18 - 0.21 cm-3. 



(n(H°)) ~0.17 cm~^ for the sightline towards the downwind cloud surface. 

Fig. 2] shows the effects of cloud opacity. The ionization levels of H, O, and 
N, which are coupled by charge exchange, decrease by >20% from the cloud 
exterior to the solar location. In contrast. He and Ne ionizations, which 
require photons more energetic by >50%, vary little. Guesstimates indicate 
that for -50% fihration of H°, converting 20% of the H from H° to H+ would 
raise the H pressure confining the paleoheliosphere at the cloud surface on the 
order of ~10%, as compared to the present-day value. The timescales over which 
these variations occur depend on the assumed cloud shape, and may be as short 
as thousands of years. 

Ionization variations that are five times larger than those we have modeled 
are observed in the very local ISM. Warm partially ionized material (WPIM, 
T > 5000 K) such as the CLIC is widespread near the Sun, and is sampled by 

^Tho small decrease in N°/H° at cloud edge results from a high N° photoionization cross- 
section. 
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Figure 4: Variation of neutral densities, due to radiative transfer effects between 
the Sun and cloud surface, for Model 2 from SF02. Shown are variations in 
neutral densities between the Sun (iV(H°)=6.5x 10^"^ cm~^) and cloud surface 
(Ar(H°)=0) for H°, He°, Ne°, 0°, and N°. At the heliopause, n(He°)~0.015 
cm^^, n{H°)^0.21 cm~'^, n(e)~0.1 cm~'^. The cloud surface is at the left, and 
the solar location is at the right of the figure. 



data on N+ and N° from Copernicus and FUSE. Values of iV(N+)/7V(N°)~0.4- 
2.0 are typical for low density sightlines, A^(H°)< 10^^ cm^^ [521121 EH]. The 
best CLIC models give n(N+)/n(N°)-0.76 at the Sun (nos. 2 and 8 in SF02). 
Consequently, as the CLIC sweeps past the Sun variations larger than 10% in 
the heliosphere may be caused by variable interstellar ionizations. The gas 
towards A Sco and HD 149499B, 37 pc away in the upwind direction, is an 
example of an ionized cloud that could engulf the Sun in the next million years. 

2.3 Solar Encounter with the CLIC 

When did the Sun enter the CLIC? With the exception of the LIC, only the 
Doppler-shifted radial components of cloud motions are observed. Hence only 
approximate estimates of the epoch that the Sun made the transition from the 
very low density Local Bubble to the CLIC are possible. We get an averaged 
value for the time of encounter by calculating the distance to the CLIC edge for 
all stars within 50 pc, and between £=170°±30° and |&| <30°, and comparing 
this distance and the HC ISM velocity towards each star. The distance to 
the cloud edge is given by A^(H°)/(n(II°)), where (n(II°))=0.17 cm~^. Nine 
stars in Figs. ^ and |21 fall in this interval. The average distance to the cloud 
edges for these stars is 2.8 pc (range 1.2-3.8 pc), and the average entry time of 
the Sun into the CLIC is 120,000 years ago (range 56,000-200,000 years ago). 
Restricting the downwind stars to those within 30 parsecs yields four stars, and 
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does not significantly change these entry time values. The star closest to the HC 
downwind direction of the LIC (^,6~183° -16°) is Ori (HD 39587, at 8.7 pc). 
The HC cloud velocity of 21.6 km s~^ towards this star [SU] indicates an entry 
date into the CLIC of ~56,000 years ago. We conclude that for no gaps in the 
ISM distribution within the CLIC, so that (n(II°))~0.17 cm"^ is vahd, then the 
Sun entered the CLIC within the past ~ 130, 000 ± 70, 000 years, and possibly 
within the past 56,000 years (consistent with our earlier estimates, |19l 124) '). 

If the velocity dispersion of the CLIC clouds represents turbulence, then the 
assumption of no gaps implies incompressible turbulence. A magnetized par- 
tially ionized tenuous plasma will have incompressible turbulence if the mag- 
netic pressure, Pb, is on the order of, or greater than, the thermal pressure, 
Pth- Magnetic field strengths of B ~3 fiG are required in the LIC to balance 
thermal pressures of Pth~2500 cm~'^ K, as expected from the best photoioniza- 
tion models with n(e)~0.1 cm'^, H+~0.1 cm'^, n(H°)~0.2 cm'^, and T=6300 
K. Weaker fields would indicate compressible turbulence and possibly gaps be- 
tween the cloudlets. In such a case, the entry of the Sun into the CLIC may 
have occurred earlier than stated above. 

The LIC is observed towards both a CMa, at 2.7 pc, and e CMa (towards 
£-230°, b — 10°, '35",T7"). The LIC iV(H°) towards a CMa has been used to set 
constraints on the entry date of the Sun into the LIC |17[ 1431 . Although the 
LIC velocity is well known, the cloud structure is uncertain, and we can only 
state that the Sun first encountered the LIC within the past 40,000 years, and 
possibly within the past 3,000 years, based on the earlier discussions. 

3 The Sun in the Local Bubble Cavity 

Prior to crossing paths with the CLIC, the Sun traveled through the Local 
Bubble (LB) for several million years The path of the Sun through the 

Local Bubble is reconstructed below from reddening data f tl3.1|l . The LB interior 
formed one type of paleoLISM, and the limits on the plasma properties of the 
Local Bubble interior are discussed in i)3.2l 

3.1 Local Bubble 

The exact dimensions and structure of the Local Bubble depend on the ISM 
component that is sampled; we use the optical reddening properties of interstel- 
lar dust grains here, based on photometric data in the Hipparcos catalog |48| . 
and find a radius of ~ 60-100 pc. Earlier maps of the Local Bubble ISM dis- 
tribution based on reddening data that traces interstellar dust pTHl ITfl 1^ IM] . 
dust and magnetic fields traced by polarization data ^UllSli and gas traced 
by H° and Na° [JS] I45L 1621 yield similar conclusions about the bubble di- 
mensions, although details of the cavity topology depend on data sensitivity to 
low column density ISM. Fig. [S] shows the motion of the Sun through the Local 
Bubble, for a plane passing through a meridian, perpendicular to the galactic 
plane, and aligned with an axis coinciding with the solar apex motion, extend- 
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ing from ^=40° to 220°. Stars with longitudes within ±25° of this meridian 
shce are plotted. Cleaned and averaged photometric and astrometric data for 
O, B, and A stars in the Hipparcos catalog are used, with variable stars and 
stars with poorly defined spectral types omitted [^Hl ■ 

In tenuous ( "intercloud" ) material, color excess, E{B — V), and 7V(H°) are 
related by logN{R)/EiB - y)>21.70 (dex cm'^ mag-\ [J). The CLIC red- 
dening is E{B — T^)<0.002 mag, and would not appear on this figure. The large 
circles show dust clouds from Dutra and Bica (2002). It is evident that for a 
solar LSR velocity of 13 — 20 pc Myrs~^, the Sun has been within the very 
low density Local Bubble for over 3 Myrs. Although CLIC column densities 
(< 10^^ cm~^) are not traced by typical E{B — V) data, ultraviolet absorption 
lines do not show any neutral ISM in the anti-apex direction, ^~220°, at the 
N(H.°)r^ 10" cm~2 igygi between ~5 and -100 pc (ED, SF06). 

3.2 Local Bubble Interior 

We now discuss the properties of the paleolism when the Sun was in the very 
low density LB interior. The nature and origin of the Local Bubble is the sub- 
ject of ongoing debate. We expect that this volume is presently filled with a 
very low density ionized plasma that provides pressure support for the cavity. 
Observations of the diffuse soft X-ray background (SXRB) led to models of the 
Local Bubble filled with hot, high pressure gas with a temperature of ~ 10^ K 
and pressure of P/fce ^ 10^ cm~'^ K {k^ is the Boltzman constant), although 
ROSAT has shown that substantial amounts of the SXRB at high latitude arises 
beyond the boundaries of the Local Bubble An effective plasma temper- 

ature near 10® K is indicated for a plasma in coUisional ionization equilibrium 
(CIE). Disregarding the problems with models of the emission spectrum, under 
the assumption that there is a CIE hot plasma filling the LB, the implied density 
is roughly 5 x 10^'^ cm"-^ and pressure of ^ P/k^ — 10"* cm^'^K. 

There are several implications of the SXRB attributed to the Local Bub- 
ble. First, the intensity of the background implies a thermal pressure that is 
substantially higher than that of the CLIC. This mismatch may be fixed by a 
CLIC magnetic field of B 4 ^G. Second, a potentially substantial contribution 
(<30% in most directions) to the SXRB may arise from charge transfer between 
solar wind ions such as O^^ and neutral H or He, either in the interstellar wind 
or geocorona [21 El ■ This does not affect the photoionization calculations de- 
scribed above directly, however, even if the heliospheric and geocoronal soft 
X-ray emission is at the upper end of current estimates. That is because the 
heliospheric and geocoronal soft X-ray emission is substantially harder than the 
emission that is directly responsible for LIC ionization (~ 13.6 — 54.4 eV). There 
is an indirect effect, though, insofar as the heliospheric X-ray emission affects 
our understanding of the source of the SXRB, which in turn has implications 
for the lower energy emission and the cloud interface. 

This low density hot LB plasma constituted the solar galactic environment 
for several million years. The LB plasma properties have probably been con- 
stant since then, because cooling times for low density plasmas are over «100 
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Myrs. Nevertheless, the LB properties are not constant. Interstellar shocks 
move through space at velocities over 100 km s~^ (~100 pc Myrs^^), quickly 
traversing low density ISM. Also, 25% of the ISM mass is contained in clouds 
traveling faster than 10 km s~^ through the LSR (based on H° 21-cm data in 
Pni). Therefore, over the past several Myrs some dilution of the plasma must 
have occurred, particularly near the LB boundaries. The CLIC, with Vlsr ~18 
km s~^ may be an example. 

The paleoheliosphere formed by the interactions between the very low den- 
sity plasma of the Local Bubble void and the solar wind has been modeled by 
Mueller et al. (2006). The resulting heliosphere dimensions are similar to the 
present heliosphere, but the heliosheath is thicker and pickup ions and other 
derivatives of interstellar neutrals are absent. These changes have opposite ef- 
fects on the modulation of galactic cosmic ray fluxes at the Earth. 

3.3 Interface between Local Bubble Plasma and CLIC 

One additional consequence of the presence of the Local Bubble plasma is that 
a transition region must exist between the tepid CLIC gas and the hot gas. 
This transition ISM constituted the paleoLISM for a very brief period of time. 
Models of the interface as a thermally conductive interface indicate a thin region 
of intermediate temperature gas with the temperature changing by an order of 
magnitude over a distance of ~ 1000 AU and more than doubling in the space of 
100 AU ((56|,SF02). The temperature profile of such a region is very steep near 
the cloud and flattens farther out as the temperature approaches that of the hot 
gas. The density proflle is similarly shaped, though inverted because the thermal 
pressure remains nearly constant in the interface. The density increases into the 
cloud by an order of magnitude over a distance of roughly 1000 AU. Over this 
same distance the ionization is decreasing from nearly complete ionization of H 
to the ~ 30% of the cloud interior. 

The cloud gas is evaporated by the influx of heat and flows off the cloud, 
accelerated by a gentle pressure gradient created by the thermal conduction. If 
such a proflle exists at the edge of the CLIC gas, the Sun would have traversed 
from gas at T ^ 10^ K to 7000 K in a period of about 500 years. This would 
result in a sudden change in the heliospheric boundary conditions that could be 
very disruptive and result in conditions far from any equilibrium conflguration. 
An alternative possibility for the boundary is that of a turbulent mixing layer 
. This would be the case if there is a substantial velocity difference between 
the CLIC gas and that of the hot gas. Even relative velocities on the order of the 
sound speed in the hot gas or less could be very disruptive to the cooler clouds 
and could lead to an interface in which the cool gas is being entrained in the 
hotter gas, mixed and then cooled. Such an interface is similar in some ways 
to the evaporative boundary described above, but depends on hydrodynamic 
instabilities to create the mixing of the cloud gas into the hot bubble gas. The 
crossing of this type of interface could also be very disruptive of the heliosphere 
with the likelihood of small scale condensations and velocity fluctuations as well 
as sudden variations in ionization in the matter incident on the Solar System. 
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4 Discussion 



Simple assumptions are necessary to estimate the eras of the solar transition 
between galactic environments, although we know that in general the ISM con- 
tains a wide range of cloud types and shows coherent structures ranging in size 
from <1 pc (e.g. the LIC) to >100 pc (e.g. Loop I). Another limitation, fol- 
lowing from the spectral analysis techniques used to acquire data, is that the 
identification of an interstellar "cloud" rests on a median velocity determined 
from the Doppler spread of atomic velocities, which is highly sensitive to the 
resolution of the instrument conducting the observations |63) . 

Over the ordinary long sightlines of the ISM, e.g. >100 pc, similar simple 
assumptions are used for clouds that are obviously blended in velocity. The 
more confined CLIC gas within ~10 pc allows higher levels of accuracy. How- 
evever, although the electron density diagnostics Mg° /Mg+ and C"''*/C+show 
consistent values of n(e)~0.1 cm~^ towards several stars sampling the CLIC, the 
neutral densities must be reconstructed from radiative transfer models because 
C*^ fine-structure data are unavailable j31j. Model results indicate consistently 
that the CLIC does not fill the sightline towards any nearby star, with filling 
factors of ~0.40 in the galactic center hemisphere and ^0.26 in the anti-center 
hemisphere for n(II°)^0.2 cm~^ j^l]. Except for observations of Si++ towards 
local stars and pickup ion Ne, radiative transfer models successfully repro- 
duce densities of ISM towards e CMa and within the solar system. Therefore, 
except for possible hidden dense small clumps of gas, we expect that a cloud 
density of n(II°)'^0.2 cm~'^ is a reasonable value. Less certain is whether there 
are gaps in the distribution of the nearest ISM, as would be expected for fila- 
mentary LIC gas |I7j. 

Nevertheless, our understanding of the nearest ISM is not complete, and in 
particular the origin of the observed macroturbulcncc indicated by the velocity 
data |22) and the possibility of denser material in the upwind direction |2L)| indi- 
cate much of the local ISM remains a mystery. The LIC is our best understood 
cloud, yet distinctly different estimates of the date of the Sun's entry in the 
LIC are found depending on which subset of data are used to define an assumed 
filamentary LIC shape (resulting in entry eras of less than 10,000 years ago 
and up to ^40,000 years ago, T71,13'). Our analysis here of the paleoLISM pro- 
vides a starting point for future studies based on improved data and a deeper 
understanding of turbulence and small-scale structure in the CLIC. 

The contribution to the soft X-ray background (SXRB) radiation field at 
the Sun due to charge exchange between the solar-wind and interstellar neu- 
trals was discussed in ^'6.'2\ The primary effect of this emission will be on the 
properties of the cloud interface, while the overall soft X-ray flux will continue 
to be dominated by the strong emission from the Loop I supernova remnant 
and attenuated extragalactic sources. 
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Table 1: HC and LSR Velocities of the Sun, CLIC, and LIC 







Sun 
Motion 


Upwind 
CLIC 


Upwind 
LIC 


HC: 

V (km 

£,b 






-28.1±4.6 
12.4°, 11.6° 


-26.3±0.4 
3.3°, 15.9° 


LSRstd : 

V (km s'^ 




19.5 


-19.4 


-20.7 


e,b 




56°, 23° 


331.0°,-5.1° 


317.8°, -0.5° 


LSRnip: 
V (km 




13.4 


-17.0 


-15.7 


e,b 




27.7°, 32.4° 


2.3°, -5.1° 


346.0°, 0.1° 



Note: The Standard solar motion, LSRstd, corresponds to a velocity of 19.5 
km s~^ towards £=56°, 5=23°. The Hipparcos solar motion, LSRnip, corre- 
sponds to a velocity of 13.4 km s"^ towards £=27.7°, 6=32.4° ,6|. 

5 Conclusions 

The primary conclusion of this paper is that, over the past several million years, 
both the galactic environment of the Sun and the heliosphere have been signif- 
icantly different than they are today. Observational data combined with the- 
oretical studies can be used to reconstruct the three-dimensional distribution 
of nearby ISM, and predict the times the Sun transitioned between different 
environments. If we assume a continuously distributed local ISM, within the 
past ~ 130, 000±70,000 years, and possibly as recent as ~56,000 years ago, the 
Sun entered low density partially ionized ISM flowing away from the direction of 
the Scorpius-Centaurus Association. Sometime within the past ~40,000 years 
the Sun entered the cloud now surrounding the solar system, the LIC. These 
estimates rely on topologically simple models of the cluster of local interstel- 
lar clouds (CLIC) flowing past the Sun; more elaborate models are discussed 
elsewhere ( \n\ 1261 SF06). As the Sun moves through this complex of local 
interstellar clouds, the boundary conditions of the heliosphere should change by 
substantial amounts due to changes in cloud temperature, velocity, and opacity- 
driven variations in the ionization of the surrounding ISM. Prior to that, the 
Sun was in the low density plasma of the Local Bubble cavity. Between the Lo- 
cal Bubble cavity and the CLIC, the Sun briefly (^500 years) passed through 
an interface region of some type. 

These estimates of the entry date of the Sun into the CLIC and LIC are 
based on current data and models. Future very-high resolution observations in 
the 1000-3000 A spectral interval, for a spatially dense sample of nearby stars, 
are required to reconstruct the distribution, kinematics, and properties of the 
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CLIC, and reduce the uncertainties in these estimates. When that happens, we 
can anticipate that the early work of Hans Fahr wiU have yielded grand results, 
as we finally understand the close relationship between the paleoheliosphere and 
paleoLISM. 
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Figure 5: The distribution of ISM is shown for non-variable stars within ±25° 
of a meridian shoe perpendicular to the galactic plane and aligned with the 
solar apex motion through the LSR. The solar apex motion is plotted as an 
arrow. The X-axis points towards a galactic longitude of -^=40°, and the Z- 
axis points towards the North Galactic Pole. The filled symbols show sightlines 
with E{B -V)>QM mag (or logA^(H)>20.48, dex cm'^), and the x's show 
sightlines with E{B -V)<Qm mag (or log7V(H)<19.93, dex cm-^). Open 
squares have E{B — y)=0.02-0.06 mag. Large open circles display the positions 
of nearby dust clouds (see text). These reddening values are determined from 
photometric data in the Hipparcos catalog j^Hj . Data points are smoothed over 
±13°, and incorporate all stars with overlapping distances (once uncertainties 
are included). The asterisks indicate the positions of the stars e CMa and /3 
CMa, which are located in a well-known spatial region with low interstellar 
densities , but which are not included in the plotted star sample because the 
Hipparcos catalog flags them as variable. 



19 



